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Abstract 

We present in this work a study of large-p-r charmonium production in hadronic col- 
lisions. We work in the framework of the factorization model of Bodwin Braaten and 
Lepage, thereby including the color octet production mechanism, and extract the values 
of the necessary nonperturbative parameters from a comparison with the most recent data 
from the Fermilab 1.8 TeV pp hadron collider. We extend the calculation to 630 GeV, and 
compare the results with data published by the UA1 Collaboration. The global agreement 
is satisfactory, indicating that the largest components of the production mechanisms for 
charmonium production at high px have been isolated. 
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1. The production of heavy quarkonium states in high energy processes has recently at- 
tracted a lot of theoretical and experimental interest. The successful operation of vertex de- 
tectors in hadronic colliders has allowed to disentangle the genuine charmonium yields from 
the large background due to production and decay of b quarks |J. The ability to detect the 
soft photons from Xc decays has recently allowed the independent measurement of the Xc con- 
tribution to the Jfip rate || [3], |]]. On the theoretical side, charmonium production provides 
quite stringent tests of our understanding of QCD on the very border between perturbative and 
nonperturbative domains. The detailed measurements of differential cross sections for produc- 
tion of J/ip, if>' and Xc states can be confronted with theoretical calculations. Starting from the 
very intuitive Colour Singlet Model ||, these have recently been improved with the inclusion 
of the mechanism of production of a parton with large transverse momentum, followed by the 
fragmentation into charmonium states |J. The inclusion of the fragmentation mechanism has 
brought the theoretical predictions closer to the observed prompt J/ip production rate @. 
However the very large discrepancy, by more than an order of magnitude, between the theo- 
retical predictions and the data for the case of the ip', clearly demands for new mechanisms 
dominating the production process. 

Several proposals have recently been put forward to solve this discrepancy. Among these, 
the possible existence of higher P-wave or D-wave states which decay into the if}', or of new 
metastable or hybrid charmonium states |§. In this paper we shall concentrate on a third 
proposal @, namely the contribution to the fragmentation function of colour octet states, 
which subsequently evolve nonperturbatively into the if)' plus soft light hadrons. 

More recently, new data on the measurement of the Xc/ip fraction |3|, |J confirm that a similar 
problem exists for J ftp' a not coming form Xc decays. The aim of the present paper is to make a 
thorough re-analysis of the full matter, trying to find a coherent picture which possibly accounts 
for the large J/i[) and ijj' production cross sections, including the new information available on 
the Xc production data. The general framework is provided by the analysis of Bodwin, Braaten 
and Lepage |10|| , which allows a consistent treatment of short and long distance effects. We will 
first review the main ingredients of this formalism, and then proceed to our phenomenological 
analysis. 

2. For the reader's convenience and to fix our notation, we briefly review some models 
which have been suggested in connection with charmonium production. We do not include a 
discussion of the color evaporation model, for which phenomenological reviews have appeared 
recently JTTJ] . We start instead presenting the Color Singlet Model (CSM) |J (for a recent 



review, see also fl^| ). The dominant mechanism is assumed to be the short- distance production 
of a color singlet QQ pair with the same spin and angular momentum quantum numbers of 
a given quarkonium state H. All the nonperturbative (long-distance) effects that lead to the 
formation of the bound state are factored into a single phenomenological parameter. Hence the 
cross section for the production of a state H = n 25+1 L j takes the form 

a[n 2S+1 Lj] = P nL a[QQ{n 25+ %)], (1) 

where the nonperturbative parameter P n i can then be expressed in terms of the radial wave 
2 Here and in the following we use the term "prompt" to refer to all sources excluding 6-decay contributions. 
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function or its derivatives, and calculated either within potential models (see for example 
ref.[[n|) or extracted from experimental data. 

This simple factorization fails in the calculation of the production of P states, for example 
via qq annihilation or in B mesons decays. In fact an infrared singularity appears, associated 
to a final state soft gluon, and at least a second nonperturbative parameter has to be invoked 
to absorb it, spoiling the simple minded picture of the CSM. 

A rigorous framework for treating quarkonium production and decays has been recently 
developed by Bodwin, Braaten and Lepage [JTTJ . Their so-called "factorization model" expresses 
the cross section for quarkonium production as a sum of terms each of which contains a short- 
distance perturbative factor and a long-distance nonperturbative matrix element, as 

^]=E^<0|<?n(A)|0> (2) 

F n are short- distance coefficients which can be calculated in perturbative QCD (pQCD) as 
power expansions in a s (m) (the quark mass m being large, a s (m) is expected to be small enough 
to allow the perturbative expansion). A is a scale which separates short and long distance effects. 
The A dependence of F n cancels against that of the matrix elements (0|O^(A)|0), leaving a 
cross section independent of A. The above matrix elements can be rigorously defined in Non 
Relativistic QCD (NRQCD) [IIJ. They absorb the nonperturbative features of the process and 
can either be extracted from data or calculated on the lattice. Finally, the 5 n are related to the 
dimension of the operator . 

The main difference between the factorization model and the CSM is that not all the oper- 
ators are now related to the production of color singlet QQ pairs. The factorization approach 
explicitly takes into account the complete structure of the quarkonium Fock space. Therefore 
the quarkonium H is no more assumed to be simply a QQ pair, but rather a superposition of 
states: 

\H = n 2S+1 Lj) = 0(l)\QQ(n 2S+1 Lj,l)) 

+ 0(v)\QQ(n 2S+1 (L±l)j,,8)g) 

+ 0(v 2 )\QQ(n 2S+1 Lj,8)gg) + ..., (3) 

where the labels 1 and 8 refer to the colour state of the QQ pair. Higher order components 
are suppressed by powers of v, the average velocity of the heavy quark in the quarkonium rest 
frame, v can be estimated through the relation 

v ~ a s {mv) , (4) 

which for the charmonium yields a value v 2 ~ 1/4. The CSM is recovered by taking only the 
lowest order term in eq. @. 

The production of the state H in the factorization approach can then proceed via any 
of the Fock components in eq. (0). Higher order components become important when their 
short distance coefficients F n are suppressed by fewer powers of ot s relative to lower order ones. 
Therefore the contribution of various terms to the production of H depends in general on both 
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the a s (m) expansion of F n and the v 2 expansion of the matrix elements: it is a two parameter 
problem. 

Similar expressions hold within the factorization model also for the fragmentation functions 
of a parton k into the state H, evaluated at a scale /i larger than the heavy quark mass: 

^(^) = £%^<0iO*(A)i0>. (5) 

The are, again, short distance coefficients. They can be calculated in pQCD at a scale Hq of 
the order of the quarkonium mass and then evolved to higher scales. After evolution, the cross 
section is given by the usual convolution: 

a[H] = f F'F^^Df (6) 



the F's being the parton distribution functions in the colliding hadrons and (Tij->k the kernel 
cross sections describing the inclusive parton-parton scattering. 

3. Within the factorization approach it is possible to relate the matrix elements of the 
leading operator in the v expansion to those entering the factorization formulae for quarko- 
nium decays. They can therefore be extracted by comparing the measured decay widths to 
those calculated. In the case of the operators relative to higher components of the Fock space 
expansion, no simple relation exists in general between decay and production matrix elements 
flT0| . So they should either be calculated (e.g. in lattice QCD), or be measured directly in some 
production process. We discuss here shortly the cases of interest for our study. More detailed 
expressions and observations can be found in fllQfl . 

In the case of Xc production, we have the following expressions for cross sections and frag- 
mentation functions to leading order in v 2 : 

a[ X j] = (0\Or( 3 Pj)\0) + (01^(^)10) J = 0,1, 2 (7) 

D f = ^T^ (0\Or( 3 Pj)\0) + dffSi) (0|^( 3 ^)|0) J = 0,1, 2 (8) 

The presence of the color octet matrix elements represents the natural extension of the CSM 
results, and allows the absorption of the infrared divergences which appear in the short distance 
coefficients of the color singlet part. Both terms are needed to give a consistent description of 
Xc production at this order in v 2 . 

The matrix elements of color singlet operators can be related to those entering the decay 
processes Xj=o,2 — ► 77 and Xj=o,2 — > light hadrons. The color octet production matrix element, 
however, cannot be related to the corresponding decay one [13]. We have to resort to a pro- 
duction process to measure it. In fl4| it was suggested to use the Xc production in B decays. 
The results have been reported in the literature in terms of the nonperturbative parameters if i 
and H$ for the color singlet and color octet parts respectively. They are related to the NRQCD 
matrix element as follows: 

i (o\or( 3 Pj )\o) 

m 4 2 J + 1 



Hi = (9) 



1 (010^(^)10) 
Hs = m 2 2J+1 (10) 
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In ref.nTJJ Hi was obtained fitting the T(x — > light hadrons), whereas in ref. M H' 8 was extracted 



from the CLEO measurement of BR(B — > Xj=i,2 + X) |15| . The values found are: 

H l « 15 MeV (11) 
F 8 = 1.4 ±0.6 MeV (12) 

Let us now consider 3 Si states, i.e. J/ip and t/>'. We will collectively indicate these as ip. The 
cross section and fragmentation function for producing a ip to leading order in v 2 are simply 
given by the CSM results: 



F 1 ( 3 S 1 ) 



<#] = ^-^(0|Of( 3 ^)|0) (13) 



ni 



D* = d^SMOfeS^O) (14) 

The matrix element appearing in the above equations can be shown to be related to the standard 
nonrelativistic wave function as follows [ID : 



(oiofCsm^^R^ 2 (15) 



They can therefore be extracted from the measurement of the leptonic decay width of the S 



states, or can be calculated within potential models |[L3j| . 

While no color octet contribution appears at order v 2 , it has recently been argued by Braaten 
and Fleming || that yet higher order terms can however be significantly enhanced since their 
short distance coefficient appears at lower orders in the a s expansion. An example of this is 
gluon fragmentation to if). To leading order in v 2 the fragmentation proceeds through a color 
singlet 3 Si state, and starts to order a s 3 : 

D^ = a s 3 D 1 (0\Ot( 3 S 1 )\0) (16) 

This is because production of the cc pair in a color singlet state requires emission of two 
perturbative gluons. On the other hand, the fragmentation process where the gluon goes into 
a color octet state, although suppressed by v 4 , starts at order a s : 

D+ = aM*\Oi{*Si)\*) (17) 

It can therefore be numerically relevant when compared to flip]). No decay process is known 
which is dominated by the color octet component, and therefore it is not possible to extract 
the relative matrix elements from decay widths. One could get a crude estimate of their values 
by rescaling the color singlet matrix elements by the appropriate powers of v. We prefer here 
to take their value as a free parameter, to be fitted to the Tevatron production data. We will 
verify at the end that the results are consistent with the expected v A suppression. 

4. We now present a comparison between experimental results and the calculations illus- 
trated above. We first concentrate on results from the Tevatron collider, relative to high-py 
production of J/ip, i>' and \c states. Since we will present results for prompt production, we will 
only make use of the CDF data, for which the 6-decay background has been removed H. It is 
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Figure 1: Inclusive prompt pr distribution. CDF data versus theory. The contri- 
bution from different sources is shown. 



important to stress, nevertheless, that there is perfect agreement between the CDF and the DO 
data when all sources of J/ip are included . We will extract the values of the nonperturbative 
parameters H%, (0g ( 3 Si)) and (Of ( 3 Si)) from fits to the experimental data. We will then use 
these values to "predict" the inclusive Jfip pt distributions at the energy of 630 GeV, where 
data are available from the measurements performed by the UAl experiment [T7I . 



For the calculation of large pt charmonium production in pp collisions at the Tevatron 
energy (y/s = 1800 GeV) we include the following contributions: 

1. Direct production of charmonium states. The matrix elements were calculated in Ref. ||. 
As previously noted in literature J7|, this contribution is very small compared to the 
fragmentation one. 

2. Production via fragmentation of gluons and charm quarks. These contributions were 
considered in ref. J/J, where it was shown that they greatly enhance the cross sections 
with respect to the direct terms. We use the fragmentation functions of gluon to ip B, 
charm to ip |L8| and gluon to x Ell- 



3. Production via fragmentation into color octet states. 
We will also show separately the contributions to J/ip and Xc production coming from the decay 



5 




4 6 8 10 12 14 

Pt f (GeV) 



Figure 2: Inclusive ip pt distribution. Upper curves and data points correspond to 
prompt ^'s, after subtraction of the Xc contribution. Lower ones correspond to the b 
decay contribution. CDF data versus theory. 



of b quarks, evaluated at the next-to-leading order [Z(| using a choice of renormalization and 
factorization scales which provides the best fit to the Tevatron data [21 . 



All the charmonium cross sections are evaluated at leading order with the MRSA [22] 
parton distribution set. The renormalization/factorization scale is set at fi = \Jpt 2 + M|- We 
have checked that using for fi the of the fragmenting parton produces differences which are 
typically of order 10-20%, therefore definitely less than the other uncertainties involved. 

Fig. [I] shows the comparison between theory and CDF data [|J for prompt ip' production. 
In this and in the following figures we have not shown the band due to theoretical uncertainties 
due to, for example, the choice of the renormalization, factorization and fragmentation scales 
j7|. This is because these uncertainties mostly affect the overall normalization of the curves, 
and not their shape. As a consequence, their effect would be hidden by a rescaling of the fitted 
value of the nonperturbative parameters. 

As was already shown in the work of Braaten and Fleming || , the theoretical curve agrees 
well with the shape of the data. The old theoretical prediction from pure color singlet frag- 
mentation was known to fall a factor of 30 below the CDF data, as shown in the figure. The 
addition of the color octet mechanism reconciles theory and data. The value we extract for 
(Of ( 3 Si)) from a best x 2 fit is 4.3 x 10~ 3 GeV 3 , close to what derived in @. 
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Figure 3: Inclusive px distribution of ^'s from Xc production and decay. Upper 
curves and data points correspond to the prompt component. Lower ones correspond 
to the b decay contribution. CDF data versus theory. 

Fig. |2] shows the inclusive px distributions of J/ip not coming from % c decays. Both prompt 
and 6-decay contributions are shown separately. The value we extract for (Of ( 3 «Si)) from a best 
X 2 fit to the prompt data is 15 x 1CT 3 GeV 3 . Without inclusion of the color octet components, 
the disagreement between theory and data would be similar that that previously noted for i/j n s, 
namely a factor of the order of 30. 

The ratio (Of ( 3 Si))/ {Of ( 3 Si)) is about 3, which is consistent with the ratio of the values 
of the color singlet wave functions at the origin. In other words, the values one extracts from 
the two independent sets of data are not unnatural within the color octet scheme. An intrinsic 
uncertainty in the comparison of these numbers comes from the ambiguity present in the choice 
of the mass scales. For example, using or M^> as opposed to 2m c in the coefficient function, 
can lead to variations up to a factor of 2 in the fit results, and in their ratios. As for the 
absolute value of the matrix elements, these are consistent with a suppression of the order of 
t> 4 ~ 0.06 relative to the color singlet ones. 

A precise prediction of the color octet mechanism, however, is that for sufficiently large 
Pt the ratio of the J/ip and ip' cross sections should be a constant. Current data do not fully 
support this expectation. It is in fact possible to obtain a good fit to the data shown in Figs. [L] 
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Figure 4: Inclusive px distribution of t/^'s at 630 GeV. All sources of ip production are 
here included. UA1 data versus theory. The parameters of the theoretical calculation 
take the values fitted on the Tevatron data. 



and E3, in the common range 4 < p? < 15 GeV, using the following parametrizations j33f : 



Mm 

dp T 
d<7(V0 



dpi 



1773exp(-p T /1.65) 
= 384exp(-p T /1.79) 



(18) 
(19) 



These fits predict a ratio which is rising with p?, varying from 0.2 to 0.4 (after removing the 
BR's) in the px range currently accessible. Taking into account the experimental uncertainties, 
this is not inconsistent with a constant ratio. If however future improved statistics should 
confirm this trend, this would be a clear indication that yet more mechanisms, such as multiple 
decays of higher charmonium resonances, are at work. 

To conclude the survey of charmonium production at the Tevatron, we present in Fig. |3] the 
Pt distribution of J/0's coming from \c decays. The theoretical curves include the effect of pt 
smearing due to the Xc Jfa decays. Both theory and data use the recent determination of 
BR{b — > Xcj + X) from CLEO [pH to extract the b — > Xc — > ip 7 contribution. 

The best x 2 n l 1° l ne prompt data gives a value of H%=3.6 MeV. The shapes of theory and 
data are consistent with each other, although the agreement is not as good as in the case of 
ip'. The value of H' 8 is larger by a factor of 2 relative to that measured by CLEO using b — > Xc 
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decays, Eq. [T2[ It should be kept in mind that the values extracted from the fits to the hadron 
collider data are directly sensitive to the perturbative K factors due to higher order corrections 
to the hard process matrix element, and to the fragmentation functions. As a reference, the 
NLO K factor for the production of a large pr gluon was evaluated to be approximately 1.5 in 
the work of Cacciari and Greco, Ref. 0. 

Having fixed the values of the nonperturbative parameters using the Tevatron data, it is 
possible to use them to make predictions for different energies and different beam types. We 
consider here data published by the UA1 Collaboration fI7 |, relative to interactions at the 630 
GeV CERN SppS Collider. We expect that at this energy and at the pt values measured by 
UAl the production mechanisms should be exactly the same as those active at the Tevatron. 
This would be true even in presence of additional processes, such as for example production 
and decay of higher charmonium resonances. 

We present these data in Fig. |], together with the theoretical predictions. UAl measured 
the J/ijj pt spectrum inclusive of all contributions from b- and Xc-decays. The contribution 
of ip' decays, once convoluted with BR(ip' — > J/ip) and with the decay spectrum, amounts to 
much less than 10% of the total, and was neglected here. The comparison shows that theory 
predicts now approximately a factor of two more J/0's than are observed. In view of what was 
said few lines above, this is contrary to our expectations. We only see one possible explanation 
for this discrepancy, leaving out experimental systematics. Namely the significant difference 
in perturbative K factors at the two energies. It has been observed since a long time that b 
production cross sections at 1.8 TeV are systematically higher than theory. Agreement with 
NLO QCD can be found only by choosing extreme values of the renormalization scale, or 
choosing values of a s larger than the input parton distribution sets prescribe, in addition to 
fixing the b mass to the relatively low value of 4.5 GeV|21|. By choosing as input parameters 
for the theoretical evaluation of the b cross section at 630 GeV the same values that fit the 
normalization of the Tevatron data, one finds a result which is approximately 30-40% higher 
than the UAl data A justification for such a discrepancy can be found in the study 

of small-x effects in heavy quark production at high energy |p5| . It is expected that such 
effects should be larger for production of c quarks, although no detailed estimate exists. If this 
were indeed the case, however, the relative discrepancy of a factor of 2 between charmonium 
production at the 1800 and at 630 GeV, as found without inclusion of small-x effects, would be 
consistent with the similar discrepancy by a factor of 30-40% found in the case of b production. 

5. We considered in this paper the large pt production of J/0, if)' and Xc states via gluon 
fragmentation into the leading color singlet and color octet components of their wave function. 
The result of a previous study by Braaten and Fleming of ip' production extends to the case of 
the J/0, showing that these effects can explain the unexpectedly large rate of prompt J/0 and 
the small Xc/ip production ratio observed at the Tevatron. The values of the nonperturbative 
parameters needed to parametrize these production processes turn out to be consistent with 
what naively expected. 

The extension of these calculations to the case of inclusive J/0 production at 630 GeV results 
in rates which are approximately a factor of 2 larger than measured by UAl. We attribute this 
discrepancy to a larger K factor at the higher energy, due to more important small- x effects. 
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Current data from the Tevatron and the residual theoretical uncertainties cannot exclude the 
presence of yet additional production mechanisms, such as the production and decay of higher 
resonances. Dominance of the production via fragmentation into the color octet component of 
the 3 Si states strictly predicts ip'/ J/xjj to be a constant, at least for pt 3> M^. Current data do 
not support this conclusions, although the statistical uncertainty is still large. 

This formalism cannot be directly applied to the calculation of total cross sections, or to 
the region pt < M$. This is because in this region the fragmentation approximation is not 
justified. The effect of color octet production, however, can be calculated including the full set 
of relevant Feynman diagrams. After this work was completed, we received a paper by Cho and 
Leibovich in which this calculation has been performed, and applied to charmonium and 
small-pr T production at the Tevatron. The results of their work are consistent with ours over 
most of the pt range covered by the charmonium data. Their value of H% is smaller than what 
we find, presumably because of the the absence in their calculation of the negative color singlet 
contribution to \ c production |Tj|. Their value of (Of ( 3 5i)) is larger than ours, presumably 
because their calculation correctly incorporates the small pt decrease in rate relative to the 
fragmentation approximation. The values of the nonperturbative parameters extracted from 
the fits to the Tevatron data can then be used, in association to the matrix elements evaluated 
25) , to perform more precise predictions of total cross sections at fixed target energies, where 



in 



a large amount of data is available. 
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